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ABSTRACT 



Aims. This work reports radiative transition rates and electron impact excitation collision strengths for levels of the 
3s 2 3p, 3s3p 2 , 3s 2 4s, and 3s 2 3d configurations of Sill. 

Methods. The radiative data were computed using the Thomas-Fermi-Dirac-Amaldi central potential, but with the 
modifications introduced by Bautista (2008) that account for the effects of electron-electron interactions. We also 
introduce new schemes for the optimization of the variational parameters of the potential. Additional calculations 
were carried out with the Relativistic Hartree-Fock and the multiconfiguration Dirac-Fock methods. Collision strengths 
in LS-coupling were calculated in the close coupling approximation with the R-matrix method. Then, fine structure 
collision strengths were obtained by means of the intermediate-coupling frame transformation (ICFT) method which 
accounts for spin-orbit coupling effects. 

Results. We present extensive comparisons between the results of different approximations and with the most recent 
calculations and experiment available in the literature. From these comparisons we derive a recommended set of gf- 
values and radiative transition rates with their corresponding estimated uncertainties. We also study the effects of 
different approximations in the representation of the target ion on the electron-impact collision strengths. Our most 
accurate set of collision strengths were integrated over a Maxwellian distribution of electron energies and the resulting 
effective collision strengths are given for a wide range of temperatures. Our results present significant differences from 
recent calculations with the B-spline non-orthogonal R-matrix method. We discuss the sources of the differences. 

Key words, atomic data - atomic processes - line: formation - stars: Eta Carinae - Active Galactic Nuclei 



1. Introduction 

Singly ionized silicon (Si n) is prominent in ultra-violet 
(UV) and optical spectra of various astrophysical plasmas. 
In terms of absorption spectra, in the spectral range long- 
ward of 912 A Sin has 8 absorption lines complexes con- 
nected to the ground term 3s 2 3p 2 P°. In photoionized plas- 
mas with electron temperatures (T e ) of the order of 10 4 K 
and electron densities (n e ) < 10 4 cm -3 the relative optical 
depths of troughs from the two ground state levels depend 
on n e . Hence, this dependence can be used as diagnostic of 
n e (e.g. Dunn et al. 2009). On this regard, the lines centered 
at 1814 A(3s 2 3p 2 P° - 3s3p 2 2 D) are particularly convenient 
because of their unusually small oscillator strengths, such 
that these, among all other Si n lines, are the most likely to 
be in the linear part of the curve of growth and their col- 
umn densities can be accurately measured. Unfortunately, 
though, the determination of accurate and reliable oscilla- 
tor strengths for these transitions is particularly difficult 
and has been the subject of much theoretical and experi- 
mental efforts. 



* Tables 6 and 7 containing the present gf-values, A-values, 
and effective collision strengths are only available in elec- 
tronic form at the CDS via anonymous ftp to cdsarc.u- 
strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi- 
bin/qcat?J/A+A/. 



Si II is also prominent in emission spectra of various 
kinds of objects. In the upper chromosphere and lower 
transition region in the Sun and late-type stars, line ra- 
tios among the Sin 1814 A multiplet and the intercombi- 
nation (3s 2 3p 2 P° - 3s3p 2 4 P) multiplet near 2335 A are po- 
tentially useful density diagnostics. However, until recently 
the best electron impact collision strengths of Dufton & 
Kingston (1991) led to predicted line ratios that disagreed 
with observations (Judge et al. 1991). Similarly, theoretical 
models of emission spectra of Broad Line Regions of Active 
Galactic Nuclei fail to reproduce the observed intensities 
by factors of a few; a problem that was named by Baldwin 
et al. (1996) the "Sin disaster". 

Despite considerable theoretical work on oscillator 
strengths there is still considerable spread in the results, 
particularly for those transitions of most astronomical in- 
terest, such as those of the 1814 A complex. This is because 
the upper 2 D term of these transitions is made of a mix- 
ture of the 3s3p 2 and 3s 2 3d configurations, which produce 
strong cancellation in the oscillator strengths and makes 
the f-values very difficult to compute. Also for these tran- 
sitions, there has been much spread in experimental gf- 
values. Measurements based on the electron beam phase 
shift method (Savage & Lawrence 1966; Curtis & Smith 
1974) could be inaccurate for very long lifetimes. Absolute 
emission measurements on an arc (Hofmann 1969) are of- 
ten uncertain owing to the difficult control and calibration 
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of the instrument. Determinations based in comparisons of 
equivalent widths of the weak 1810 A lines to stronger lines 
in astronomical spectra (Shull, Snow & York 1981; Van 
Buren 1986) are unreliable due to blends and saturation of 
the stronger lines (Jenkins 1986). The latest, and proba- 
bly most accurate determination of gf-values for these lines 
was done with the time-resolved laser-induced fluorescence 
technique Bergeson & Lawler (1993). 

The oscillator strengths for the intercombination tran- 
sitions (3s 2 3p 2 P° - 3s3p 2 4 P) have also been subject of 
controversy. These transitions arise due to spin-orbit mix- 
ing of the metastable term 3s3p 2 4 P with the 2 S, 2 P, and 
2 D terms. The difficulty then comes from the fact that the 
gf-values are sensitive to core-valence and core-core corre- 
lations, in addition to the valence-shell effects. The results 
of some of the most recent calculations (Froese Fischer et 
al. 2006) differ from experimentally measured transition 
probabilities (Calamai et al. 1993) by ~ 30%. 

But beyond providing accurate oscillator strengths by 
tailoring atomic structure representations on each type of 
transition of interest, it is important to try constructing a 
single representation that yields reasonably accurate ener- 
gies and oscillator strengths for all levels and transitions 
considered simultaneously. This is because various practi- 
cal applications will require such a general atomic repre- 
sentation for subsequent scattering calculations (e.g. elec- 
tron impact excitation, photoionization, and recombina- 
tion). One of the most widely used methods for this purpose 
is based on the Thomas-Fcrmi-Dirac-Amaldi (TFDA) cen- 
tral potential to generate optimized one-electron orbitals, 
which represent the atomic structure through configura- 
tion interaction (CI) expansions and can also be used for 
close-coupling representations of the scattering problem. 
Recently, Bautista (2008) introduced a correction to the 
TFDA potential that accounts for the effects of electron- 
electron interactions on the radial wavefunction and yield 
a considerably improved representation of the system. It is 
thus interesting to apply this new approach to get the best 
possible representation of the important Si n system. 

In the present paper, to complement and expand the 
work with the TFDA potential we also compare with the re- 
sults of the Hartree-Fock Relativistic (HFR) , and the multi- 
configuration Dirac-Fock (MCDF) methods. This general 
multiplatform approach was successfully employed in our 
previous studies of the K-shell spectra of Fe, O, Ne, Mg, Si, 
S, Ar, Ca, and Ni (e.g. Bautista et al. 2003, Garcia et al. 
2005, Palmeri et al. 2003a, 2003b, 2008a, 2008b). This has 
the advantage that allows for consistency checks and inter- 
comparison. It also helps to reveal which physical processes 
are important for any given transition, since these different 
platforms employ different approaches to, for example, rel- 
ativistic effects and orthogonal vs. non-orthogonal orbitals. 

The present paper is organized as follows: In the next 
section we describe the calculations of oscillator strengths 
and transition rates. In Section 3 we present the calcula- 
tions of collision strengths, including some new techniques 
specifically designed for the present problem. Our discus- 
sion of the results and conclusions are presented in Section 
4. 



2. Radiative Calculations 

The Brcit Pauli Hamiltonian for an A r -clectron system is 
given by 

Hbp = H nl + Hi\y + H 2 b (1) 

where H nl is the usual non-relativistic Hamiltonian, and 
-ffib and H 2 h are the one-body and two-body operators. 
The one-body relativistic operators 

N 

H lh = ]T /„(mass) + /„(d) + /„(so) (2) 

n=l 

represent the spin-orbit interaction, /„(so), and the non- 
fine structure mass- variation, /„(mass), and one-body 
Darwin, / n (d), corrections. The two-body corrections 

#2b = 9n,n{so) + 9n m (ss) + #„ TO (cSs) + g nm (d) + 

+g nm (oo) , (3) 

usually referred to as the Breit interaction, include, on the 
one hand, the fine structure terms g n m(so) (spin-other- 
orbit and mutual spin-orbit) and g nm (ss) (spin-spin); and 
on the other, the non-fine structure terms: g nm (css) (spin- 
spin contact), g nm (d) (Darwin) and g n m(oo) (orbit-orbit). 

The oscillator strengths (/-values) for dipole allowed 
transitions have equivalent forms in length and velocity 
gauges as 



fjk 



2m(E k -E j ) 



(|j|x|fc)| 2 = ^|(|j|v|. 



(4) 



The radiative rates (A-values) for electric dipole and 
quadrupole transitions are respectively given in units of s _1 
by the expressions 

A E1 (k,i) = 2.6774 x l0 9 (E k - E.f — SM^i) (5) 

9k 

A E2 (k,i) = 2.6733 x l0 3 (E k - E t f — S E2 (M) (6) 

9k 

where S(k, i) is the line strength, g k is the statistical weight 
of the upper level, and energies are in Rydberg units and 
lengths in Bohr radii. 

Similarly for magnetic dipole and quadrupole transi- 
tions, the A-values are respectively given by 

A M1 (k, i) = 3.5644 x 10 4 (E k - E t ) 3 — S M i(k, i) (7) 

9k 

A M2 {k,i) = 2.3727 x HT 2 (£ fe -£,) 5 — S M2 (M) • (8) 

9k 

Due to the strong magnetic interactions in this ion, the 
magnetic dipole line strength is assumed to take the form 

s M i(M)H(|fc|P|z)| 2 (9) 

where 

N 

P = P° + P 1 = ^{l(n) + ( j(n)} + P rc . (10) 

n=l 

P° is the usual low-order Ml operator while P rc includes 
the relativistic corrections established by Drake (1971). 

In the present work we employ three different compu- 
tational packages to study the properties of the strongly 
correlated configurations 3s 2 3p, 3s3p 2 , 3s 2 3d, and 3s 2 4s of 
Sin. 
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2.0.1. AUTOSTRUCTURE 

AUTOSTRUCTURE, an extension by Badnell (1986, 1997) of 
the atomic structure program SUPERSTRUCTURE (Eissner 
et al. 1974), computes fine-structure level energies, radia- 
tive and Auger rates in a Breit-Pauli relativistic frame- 
work. Single electron orbitals, P n i(r), are constructed by 
diagonalizing the non- relativistic Hamiltonian, iJ nr , within 
a statistical Thomas-Fermi-Dirac-Amaldi (TFDA) model 
potential V{\ n i) (Eissner & Nussbaumer 1969). The X n i 
scaling parameters are optimized variationally by minimiz- 
ing a weighted sum of the LS term energies. LS terms 
are represented by configuration-interaction (CI) wavefunc- 
tions of the type: 



(ii) 



Bautista (2008) introduced a modification to the TFDA 
model potential that accounts in part for the effects of 
electron-electron correlations on the radial wavefunctions 
by means of additional higher order terms in the potential 
whose strength is controlled by variational parameters. We 
will refer to this potential as c-TFDA. This formalism also 
proposed a new optimization technique of the variational 
parameters by minimizing the difference between predicted, 
including spin-orbit coupling and relativistic effects, and 
experimental term averaged energies. The numerical func- 
tional used for this optimization was: 



\Ef s -E\ heo \ +e 
E obs 



(12) 



with E° bs and Ef heo the observed and theoretical energies 
respectively of term i and 

e = max{0, (E core - E° ore )} 

This last term acts as a penalty to the functional whenever 
the new core energy exceeds that found from the standard 
TFDA minimization. This numerical functional was modi- 
fied again for the present work as shown in Section 2.2. 



2.0.2. HFR 

In the Hartree-Fock Relativistic code (hfr) by Cowan 
(1981), a set of orbitals are obtained for each electronic 
configuration by solving the Hartree-Fock equations for 
the spherically averaged atom. The equations are the re- 
sult of the application of the variational principle to the 
configuration average energy. Relativistic corrections are 
also included in this set of equations, i.e. the Blume- 
Watson spin-orbit, mass- variation and one-body Darwin 
terms. The Blume-Watson spin-orbit term comprises the 
part of the Breit interaction that can be reduced to a one- 
body operator. 

The multi-configuration Hamiltonian matrix is con- 
structed and diagonalizcd in the LSJir representation 
within the framework of the Slater-Condon theory. Each 
matrix element is a sum of products of Racah angular coeffi- 
cients and radial integrals (Slater and spin-orbit integrals) , 
i.e. 



The radial parameters, , can be adjusted to fit the avail- 
able experimental energy levels in a least-squares approach. 
The eigenvalues and the eigenstates obtained in this way 
(ab initio or semi-empirically) are used to compute the 
wavelength and oscillator strength for each possible transi- 
tion. 

2.0.3. MCDF 

The multi-configuration Dirac-Fock (MCDF) method con- 
siders the Dirac Hamiltonian for a N-electron atomic system 
given in a.u. by: 



H 



D 



N 

E 

i=l 



ca\ x pi - (pi - l)c 2 



Z 
r, 



(14) 



(a\H\b) 



Ea.b T a,b 



(13) 



where c is the speed of light and a and (3 are the Dirac 
matrices. The atomic state function (ASF) is given as 
an expansion over jj-coupled configuration state functions 
(CSFs), which in turn are constructed from Slater de- 
terminants built on four-component Dirac orbitals. The 
MCDF method is implemented in the computer package 
GRASP2K, described by Jonsson et al. (2007). 

2.1. Energy levels and radiative rates 

Our calculations concentrate on the 15 lowest levels, of the 
Si II system. This is 8 LS terms from configurations 3s 2 3p, 
3s3p 2 , 3s 2 4s, 3s 2 3d, and 3s 2 4p. A Grotrian diagram of the 
ions is presented in Fig. 1 which shows the levels and the 
UV optical and intercombination transitions of interest. 

2.1.1. Calculations with HFR 

In the physical model adopted in HFR calculations, we sup- 
pose that the Si + ion can be represented by three valence 
electrons surrounding a Nc-like Si 4 + ionic core with ten 
electrons occupying the ls 2 2s 2 2p 6 closed subshells. The in- 
travalence correlation is then considered by the explicit in- 
troduction, in the model, of the 52 following configurations 
: 3s23p + 3s24p + 3s25p + 3s 2 4f + 3s 2 5f + 3s3p3d + 
3s3p4d + 3s3p5d + 3s3p4s + 3s3p5s + 3s3d4p + 3s3d5p + 
3p 2 4p + 3p 2 5p + 3p 2 4f + 3p 2 5f + 3d 2 4p + 3d 2 5p + 3d 2 4f 
+ 3d 2 5f + 3p 3 (odd parity) and 3s3p 2 + 3s 2 4s + 3s 2 5s 
+ 3s 2 3d + 3s 2 4d + 3s 2 5d + 3s3p4p + 3s3p5p + 3s3p4f + 
3s3p5f + 3s3d 2 + 3s4s 2 + 3s4p 2 + 3s4d2 + 3s4s5s + 3s3d4s 
+ 3s3d5s + 3s3d4d + 3s3d5d + 3p 2 3d + 3p 2 4d + 3p 2 5d 
+ 3p 2 4s + 3p 2 5s + 3p3d4p + 3p4s4p + 3d 2 4s + 3d 2 5s 
+ 3d 2 4d + 3d 2 5d +3d 3 (even parity). Core- valence corre- 
lation is then considered by including a core-polarization 
(CPOL) potential and a correction to the dipole operator 
as described in many previous papers (see e.g. Quinet et 
al. 1999, 2002). These corrections are used with a value of 
the dipole polarizability equal to 0.16 ai 3 ,, as computed by 
Johnson et al. 1983 for Si V, and a cut-off radius equal to 
0.53 ao which corresponds to the HFR expectation value of 
< r > for the outermost core orbital, i.e. 2p. This method 
is then combined with a well-established least-squares opti- 
mization of the radial parameters in order to minimize the 
discrepancies between the Hamiltonian eigenvalues and the 
available experimental energy levels for the 3s23p, 3s 2 4p, 
3s3p 2 , 3s 2 4s, 3s 2 5s, 3s 2 3d, 3s 2 4d, 3s 2 5d and 3s3p4p config- 
urations. 



4 



M.A. Bautista et al.: Atomic data from the Iron Project LXVI 




Fig. 1. Partial Grotrian diagram of Sill showing the 22 
lowest Levels in 8 LS terms and the dipole and intercombi- 
nation line multiplets observable in the UV spectral region. 

2.1.2. Calculations with MCDF 

Here, we carried out two MCDF calculations. The first one, 
hereafter referred to as MCDF1, was focused on our radia- 
tive parameters for the intercombination transitions 3s 2 3p 
2 P° - 3s3p 2 4 P, and the second one, hereafter referred to 
as MCDF2, on those of the allowed transitions 3s 2 3p 2 P° - 
3s3p 2 2 D. 

In MCDF1, CI expansions are built including the 
valence- valence correlation through single and double ex- 
citations from 3s 2 3p J=l/2,3/2 + 3s3p 2 J=l/2,3/2,5/2 ex- 
tending the orbital active set up to 5g in 4 steps. The final 
number of CSFs generated was 2127. For 3 of the steps, 
we selected an "extended optimal level", EOL, optimiza- 
tion option (Dyall et al. 1989; Parpa et al. 1996) on the 
lowest 5 levels using the same weight for the 5 ASFs. In 
the first step, no excitation was allowed, the active set was 
limited up to 3p, and all the core and valence orbitals were 
variational. The second step consisted in extending the ex- 
pansion through single and double excitations increasing 
the active set up to 3d. Here, only the 3d orbital was vari- 
ational, while freezing all the others to their values of the 
preceding optimization step. In the third step, the expan- 
sions were further extended increasing the active set up to 
4f. All n=4 orbitals were optimized while freezing all oth- 
ers in a similar procedure as in the second step. In the final 
optimization step, only the n=5 were variational freezing 
all the other orbitals to their preceding step values. 



In MCDF2, we used the same procedure as in MCDF1 
except that the excitations were from 3s 2 3p J=l/2,3/2 + 
3s3p 2 J=3/2,5/2 and the EOL option was selected on the 
first, second, fifth and sixth levels. Here the final number 
of CSFs used in the MCDF expansions was 1832. 

2.1.3. Calculations with autostructure 

We performed calculations with various different configu- 
rations expansions, starting with models similar to those of 
previously published work and then evolving to more so- 
phisticated techniques and larger configuration expansions. 

Table 1. Configurations included in various expansions 
used for Sill 



CE1 2s :4 2p b 3s i! 3p, 2s^2p b 3s3p^, 2s^2p b 3s3p3d, 
2s 2 2p 6 3s3d 2 , 2s 2 2p 6 3p 3 , 2s 2 2p 6 3p 2 3d, 
2s 2 2p 6 3p3d 2 , 2s 2 2p 6 3d 3 , 2s 2 2p 6 3s 2 3d, 
2s 2 2p 6 3s 2 4s, 2s 2 2p 6 3s 2 4p, 2s 2 2p 6 3s 2 4d, 
2s 2 2p 6 3s 2 5s, 2s 2 2p 6 3s 2 5p, 2s 2 2p 6 3s 2 5d, 
2s 2 2p 6 3s3p4s, 2s 2 2p 6 3s3p4p, 2s 2 2p 6 3s3d4p, 
2s 2 2p 6 3s3d4d, 2s 2 2p 6 3s4p 2 , 2s 2 2p 6 3s4s5s, 
2s 2 2p 6 3p 2 4s, 2s 2 2p 6 3p3d4p, 2s 2 2p 6 3d 2 4s 



CE2 


2s ;i 2p b 3sMf, 2s^2p D 3s 2 5f 


CE3 


2s i! 2p ;> 3s^3p^, 2s^2p & 3s^3p3d, 2s z 2p & 3s ;i 3p4s 


CE4 


2s ;i 3p 4 3s ; '3p J , 2s ;i 2p 4 3s ;i 3p^3d, 2s ;i 2p 4 3s ;i 3pMs 


CE5 


2s2p b 3s 2 3p 2 , 2s2p b 3s 2 3p3d, 2s2p b 3s^3p4s 


CE6 


2s 2 2p t> 3s 2 3p4p, 2s 2 2p t> 3s 2 3p4d, 2s 2 2p 4 3s 2 3p 2 4p, 




2s 2 2p 4 3s 2 3p 2 4d, 2s 2 2p 3 3s 2 3p 4 , 2s 2 2p 3 3s 2 3p 3 3d, 




2s 2 2p 3 3s 2 3p 3 4s, 2s 2 2p 3 3s 2 3p 3 4p, 2s 2 2p 3 3s 2 3p 3 4d 



Table 1 presents the configurations included in 7 basic 
expansions studied. The various configuration expansions 
(CEs) are named CE1, CE2, CE6 and are tabulated 
in such a way that each CE includes the configurations 
of all previous CEs. In other words, CE1 includes the 24 
configurations listed at the top of Table 1, CE2 includes all 
24 configurations from CE1 plus 2s 2 2p 6 3s 2 4f, CE3 includes 
the 25 configuration from CE1 and CE2 plus 3 more, and 
CE6 includes all 48 configurations in the table. 

The rational for these configurations and the progres- 
sion of various calculations carried is as follows: 
- AST1: This calculation uses the expansion CE1 that 
keeps the 2s and 2p shells closed and promotes the three 
remaining electrons among orbitals with principal quan- 
tum number n=3, 4, and 5 and orbital angular momentum 
/=0, 1, 2. The TFDA potential was used and the orbitals 
were optimized in the standard AUTOSTRUCTURE proce- 
dure minimizing the energies of the lowest 8 LS terms. Fine 
structure coupling and relativistic corrections are intro- 
duced as perturbations after optimizing the orbitals. From 
this calculation the core energy of the system is -578.40 Ry. 
Table 1 compares the calculated energies relative to the 
ground level with experimental values from NIST (2008). 
The results are rather unsatisfactory with energy differ- 
ences scattered between -19% and 35%. Even more trou- 
blesome are the q f- values that result from this calculation, 
particularly those for the 1814A multiplet, which are con- 
siderably overestimated with respect to the experimental 
determination of Bergeson & Lawler (1993) (see Fig. 2). 
This is because it is difficult to represent 3s 2 3p and 3s3p 2 



M.A. Bautista et al.: Atomic data from the Iron Project LXVI 
Table 2. Comparison of level energies (Ry) for Sill 



5 



1/2 



r l/2 
2po 

r 3/2 
4p 

4 P 3 /2 
JP 5 /2 

2 D 5/2 



3s3p 2 
3s 2 3d 
3s 2 3d 
3s 2 4p 
3s 2 4p 
3s3p 2 
3s3p 2 



0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 

0.002618 -18.8 -17.8 -14.8 -15.0 -15.0 -15.0 -15.1 0.19 -0.12 

0.390244 -7.86 -17.4 -11.3 -11.3 -10.9 -10.9 -10.6 -7.06 -7.62 

0.391231 -7.90 -17.4 -11.4 -11.3 -10.9 -10.9 -10.6 -7.01 -7.57 

0.392828 -7.96 -17.4 -11.4 -11.3 -10.9 -10.9 -10.6 -6.92 -7.48 

0.504016 14.5 2.42 1.03 1.00 1.18 1.18 1.42 -0.01 -0.25 

0.504160 14.4 2.40 1.01 0.98 1.15 1.16 1.39 -0.02 -0.27 

0.596884 -0.53 3.31 -0.39 -0.49 -0.59 -0.54 -0.32 15.3 -4.64 

0.698626 0.69 1.64 3.23 3.25 3.46 3.46 3.60 2.73 1.78 

0.722986 35.3 -2.37 1.19 1.13 1.14 1.17 1.29 11.1 10.3 

0.723136 35.2 -2.39 1.17 1.12 1.12 1.16 1.28 11.1 10.3 

0.739870 22.0 3.40 6.77 6.83 7.03 7.07 6.81 13.9 13.7 

0.740416 22.3 3.41 6.82 6.87 7.08 7.12 6.86 14.0 13.9 

0.763660 0.75 5.74 5.52 5.49 5.68 5.71 5.74 4.20 5.01 

0.765503 0.69 5.68 5.47 5.44 5.64 5.66 5.70 4.21 5.04 



states simultaneously with orthogonal orbitals, owing to po- 
larization and orbital relaxation effects. 

- AST2: For this calculation we use the same expansion 
CE1 as in AST1 but with the use of the c-TFDA potential 
of Bautista (2008) optimized against the experimental term 
energies using the functional of Eqn. 2.0.1. This gave a core 
energy of -578.34 Ry which is slightly higher than in AST1. 
Yet, the calculated energies here agree much better with 
experiment than from AST1, with most terms within 5% 
of experiment with the exception of the two excited levels 
3s 2 3p 2 ?3/27 fr° m the ground multiplct, and 3s3p 2 4 Pj. 
Together with the changes in calculated energies there are 
also significant changes in the predicted ^/-values for the 
1814A and 153lA multiplets, as illustrated in Fig. 2. 

- AST3: This calculation is just like AST2, but we now 
add the configurations 3s 2 4f and 3s 2 5f (i.e. CE2). We find 
that the f orbitals have an important effect on increasing the 
polarizability of configurations that include d orbitals, such 
as the 3s3d 2 , and 3p 2 3d. Thus, by including nf configura- 
tions in the expansion we find different optimization for the 
core. When including these ni orbitals one has to include 
enough nd configurations in order to reach convergence in 
the solution. But, various configurations such as 3d, 3d2, 
4d, and 5d are already included since CE1. In the course 
of the calculations we checked on the contributions of 5d 
and 5f and they were found to be small. Thus, we believe 
that the calculation is nearly converged. In this calculation 
the core energy improves to -578.38 Ry and the agreement 
between calculated and experimental energies with respect 
to the ground level also improves a little. Despite the rela- 
tively minor changes in the calculated energies between this 
calculation and the previous one there are large changes in 
the ^/-values for the 1814A and 1264A multiplets, with a 
drop of ~0.8 dex for the former and a rise of ^0.9 dex for 
the latter. 

- AST4, AST5, AST6, and AST7: For these calculations 
we progressively add new configurations by opening the 2s 
and 2p orbitals, with one 2p electron promotions in AST4 
(CE3), two 2p electron promotions in AST5 (CE4), one Is 
electron promotions in AST6 (CE5) and three 2p electron 
promotions in AST7 (CE6). By this point the computa- 



tional time in optimizing all orbitals of the potential has 
become exceedingly long, going from a few minutes in the 
AST1 calculation to several hours in the AST7 calculation. 
For this reason, and with the purpose to see the effects of in- 
creasing CI we did not optimize again the orbitals for these 
calculation, but instead chose the same scaling parameters 
as in AST3. Despite the large increase in CI in these com- 
putations we find no significant change in either the cal- 
culated energies, see Table 2. The calculated core energies 
also remain constant at -578.38 Ry for AST4, AST5, and 
AST6 and -578.39 Ry for AST7. The computed /-values 
are also unaltered by CI, except for those of the 1814A mul- 
tiplet which seem to converge asymptotically. Nonetheless, 
for this and most of the other transitions the difference be- 
tween length and velocity gauges of the /-values is uncom- 
fortably large. It is clear now that increasing the amount 
of CI in the calculations beyond that of the AST7 calcula- 
tion will not improve the quality of the results. Instead, we 
need to turn our attention towards the optimization of the 
orbitals. 

- AST8. For this calculation we decided to return to the 
simplest expansion CE1, but try a different optimization 
technique. Because the /-values for the 1814A multiplet is 
so sensitive to cancellation in the mixing of the 3s3p 2 2 D 
and 3s 2 3d 2 D states it is crucial to reproduce the experi- 
mental energies of the levels as accurately as possible. Thus, 
we modified the AUTOSTRUCTURE code use an optimization 
functional of the form 

( T?obs rptheo\2 i 

^ WjjEj -Ej ) + 
r (E obs ) 2 

where the sum goes over fine structure energy levels and wt 
is a user defined weight on different levels. For our calcu- 
lation we ran this sum for the lowest 15 levels of the Si II 
system and we weighted the 3s 2 3p 2 Pg/ 2 an d 3s3p 2 2 Dj 
levels 5 times as much as the rest. This kind of optimiza- 
tion provides much better energies for all levels and nearly 
exact energies for the 3s 2 3p 2 Pg/ 2 and 3s3p 2 2 Dj levels 
(within 0.2% of experiment). This optimization also pro- 
vides a lower core energy than from previous calculations 
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Fig. 2. Evolution of the (//-values for transitions in Sin 
from our different calculations AST1, AST9 and com- 
parison with previous determinations. The square points 
depict our results in the length gauge (blue) and velocity 
gauge (red). The solid horizontal lines depict experimental 
determinations by Shull, Snow & York (1981) (magenta), 
Van Buren (1986) (blue), Schectman et al. (1988) (green), 
Bcrgeson & Lawler (1993) (black). The dashed lines depict 
theoretical results of Froese Fischer et al. (2006), (blue), 
Tayal (2007) (red), Nahar (1988) (green). 

(-578.42 Ry). Nevertheless, this optimization does not re- 
duce the difference between length and velocity forms of 
the /-values. 

An alternative optimization functional that is already 
available in AUTOSTRUCTURE is based on the differences 
between length and velocity /-values. This, however, does 
not work because it leaves the level energies to change way 
off from experiment. 

Then, we define a combined functional of the form 



F = 



E 



Wi (E? 



E 



theo\2 



) 2 + e 



obs\2 



gfZf 



(gf, 



9f? 3 ) 2 



(16) 



where gf\^ and gf^ are the length and velocity forms of 
the gf- value for the transition i — ► j. An optimization of 
this kind heavily weighted on the ^/-values for the 1260A 
multiplet yields length and velocity values in very close 
agreement with each other. These results also agree very 



Fig. 3. Cont. Fig. 2 



well with the experimental determination of Bergeson & 
Lawler (1993). Unfortunately, the obtained gf- values for 
transitions involving higher terms seem to deteriorate. This 
is to be expected, as the small expansion used here misses 
important core-valence and core-core correlations. 
- AST9: For this last calculation we use the largest ex- 
pansion CE6 and optimize the orbitals as in AST8. For 
a model expansion of this size numerical optimization be- 
comes challenging, and the code took 4 days to do so. The 
resulting energies are slightly better than in the previous 
calculation, and the core energy is predicted at -572.41 Ry. 
It is found that the /-values in the length gauge for the 
1814A multiplet remain in good agreement with the ex- 
periment of Bergeson & Lawler (1993) although the ve- 
locity /-value departed somewhat. Interestingly, for these 
transitions the present velocity /-value is lower than the 
length /-value, contrary to what was seen in previous cal- 
culations. The present difference between the two forms of 
the /-value probably owes to the need of finer optimization 
of the orbitals, which became computationally prohibitive 
for the large CI expansion used here. Another interesting 
set of /-values are those for the 1260A multiplet, which in 
the present calculation result somewhat higher and with 
greater difference between length and velocity than from 
previous calculations. The present result seems to favor the 
experimental determination of Schectman et al. (1988), 
while the results of previous calculations seemed very sta- 
ble and in good agreement with Froese Fischer et al. (2006) 
and Tayal (2007). 
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Fig. 4. Evolution of the gf- values for intercombination tran- 
sition in Si II from our different calculations AST1, AST9 
and comparison with experimental values of Calamai et al. 
(1993). The square points depict our results in the length 
gauge (blue) and velocity gauge (red). The experimental 
vales are given by the horizontal solid line, with quoted 
errors given by dashed lines. 



At this point is worth pointing out the detailed study 
of Dufton et al. (1992) on the LS f-value for the 1814 A 
transition. They found that in dealing with the cancelation 
effects from contributions from the 3s3p 2 2 D and 3s 2 3d 2 D 
states one must pay attention to the calculated energy dif- 
ference between these states. In this sense we notice that 
our prefer calculation AST8 and AST9 yield energy sep- 
arations for the 3s3p 2 2 D and 3s 2 3d 2 D states that differ 
by ~ 10 — 11% from experiment, in contrast with calcu- 
lations AST3 through AST7 that reproduce these energy 
separation significantly better. However, getting exact en- 
ergies is neither sufficient nor guarantee of resulting accu- 
rate gf-vaues. Inspection of Table 1 of Dufton et al. (1992) 
shows that essentially the same gf-value can obtained when 
the energy difference with respect to experiment was less 
than 2% and nearly 11%. By contrast, calculations that pre- 
dicted energies differing from experiment by 5 and 7% got 
gf-values that were off by roughly a factor of 3. Moreover, 
we are mostly concerned with the fine structure gf-values 
and for that we also need accurate energies for the J re- 
solved levels of the 3s 2 3p 2 P and 3s3p 2 2 D terms and these 
are best reproduced by the AST8 and AST9 calculations, 



as opposed to the AST7 calculation that underestimates 
the energy of the 3s 2 3p 2 P3/2 level by 15%. 

On the other hand, it is also known that in comput- 
ing gf-values one may not rely only on the agreement be- 
tween length and velocity values, because such an criterium 
can be misleading if the model potential is oversimplified. 
Nonetheless, we believe that the agreement between length 
and velocity results in our calculations AST8 and, partic- 
ularly AST9 is very significant. AST9 uses the largest CI 
expansion that we could handle, which the inner core or- 
bitals to account for polarization. In addition, this calcula- 
tion uses the modified TFDA potential to account for ad- 
ditional electron-electron correlation. Thus AST9 employs 
the most complex and complete potential that we can build 
at this time, then the good agreement between length and 
velocity seems meaningful. 

Now we look at the results for the intercombination 
transitions 3s 2 3p 2 P? - 3s3p 2 4 Pj'. These lines, seen in 
emission in the 2328 - 2350A spectral range, are very 
useful density diagnostics in a variety of astronomical 
sources. Though, computation of accurate transition rates, 
A-values, for these lines is difficult because the transitions 
are induced by mixing between the 4 P and even parity dou- 
blet states. Hence, the calculated A-values for the intercom- 
bination transitions depend critically on the quality of the 
representation of those even parity doublet states. Fig. 4 
shows the evolution of the calculated (//-values in length 
and velocity forms from our calculations together with the 
results of other authors and the experimental determina- 
tions of Calamai et al. (1993). Here we plot our results 
in the length and velocity gauges, although none of previ- 
ous authors present similar comparisons. It is clear, though, 
that there are considerable discrepancies between the two 
gauges, being the length form of the (//-values typically 
in better agreement with experiment. Other authors, like 
Tayal (2007) do not even quote their results in the velocity 
gauge. Indeed, the velocity form of the (//-value is typi- 
cally less stable numerically, since it depends on the first 
derivative of the radial wavefunction. Thus, we disregard 
the velocity ^/-values for subsequent analysis. The present 
(//-values from AUTOSTRUCTURE agree within 0.2 dex with 
experimental the determinations of Calamai et al. (1993), 
and the agreement between MCDF results and experiment 
is even better, within 0.1 dex. On the other hand, the HFR 
results seem as accurate as MCDF for some transitions, 
but significantly discrepant for others. The results of Tayal 
(2007) and Froese Fischer et al. (2006), using a multicon- 
figuration Hartree-Fock, are comparable in quality to our 
HFR results. The results of Dufton & Kingston (1991), 
based on the multiconfiguration Hartree-Fock method, and 
Nussbaumer (1977), from the use of the TFDA potential, 
seem to be of inferior quality. 

It seems clear that no one calculation among those per- 
formed here and those reported by other authors can pro- 
vide ultimate accuracy (//-values for all transitions simul- 
taneously. Like with the different calculations, there is sig- 
nificant scatter among the results of various experimental 
determinations. Thus, in order to provide the most reli- 
able set of (//-values possible we take the statistical average 
among all theoretical results, i.e. our present results from 
AST8 and AST9 and those of other authors. In computing 
the average we discard values that depart by more than 3a 
from the average, where a is the statistical dispersion of the 



Table 3. logfp/) values for transitions among the 15 lowest levels in Sill 



Lower level 


Upper level 


A (A) 


AST8 


AST9 


HFR 


MCDF 


Rccom. 


Error 






log(s/)i 


log(9/)„ 


l°g(9/)( 


log(s/)„ 


log(s/) 


l°g(9/)c 


log(s/)s 






3s 2 3p "p; /2 


3s3p 2 4 P 1/2 


2334.41 


-5.00 


-4.89 


-4.89 


-4.55 


-4.98 


-5.09 


-5.10 


-5.05 


0.14 


3s 2 3p 2 P° /2 


3s3p 2 4 P 1/2 


2350.17 


-5.08 


-4.66 


-4.94 


-4.86 


-5.13 


-5.11 


-5.13 


-5.18 


0.15 


3s 2 3p 2 P° /2 


3s3p 2 4 P 3/2 


2328.52 


-7.72 


-7.28 


-7.49 


-8.06 


-7.19 


-6.53 


-6.62 


-7.39 


0.24 


3s 2 3p 2 P° /2 


3s3p 2 4 P 3/2 


2344.20 


-5.32 


-4.42 


-5.27 


-4.59 


-5.15 


-5.43 


-5.44 


-5.39 


0.16 


3s3p 2 4 P 5/2 


2334.20 


-5.02 


-4.20 


-4.94 


-4.40 


-4.80 


-4.81 


-4.85 


-4.99 


0.20 


3s 2 3p 2 P° /2 


3s3p 2 2 D 3/2 


1808.00 


-2.39 


-2.33 


-2.47 


-2.90 


-2.49 


-1.99 


-2.03 


-2.31 


0.16 


3s 2 3p 2 P 3/2 


3s3p 2 2 D 3/2 


1817.45 


-3.34 


-3.34 


-3.45 


-4.28 


-3.45 


-2.82 


-2.87 


-3.14 


0.26 


3s3p 2 2 D 5/2 


1816.92 


-2.20 


-2.10 


-2.30 


-2.72 


-2.28 






-2.17 


0.22 


3s 2 3p 2 P° /2 


3s 2 4s 2 S 1/2 


1526.72 


-0.24 


-0.38 


-0.46 


-0.51 


-0.60 






-0.54 


0.10 


3s 2 3p 2 P 3/2 


3s 2 4s 2 Si /2 


1533.45 


+0.05 


-0.11 


-0.46 


-0.51 


-0.30 






-0.27 


0.09 


3s 2 3p 2 P°, 2 


3s3p 2 2 Si/ 2 


1304.37 


-2.62 


-2.07 


-0.71 


-1.01 


-0.64 






-0.74 


0.09 


3s 2 3p 2 P° /2 


3s3p 2 2 Si/ 2 


1309.27 


-3.12 


-2.83 


-0.47 


-0.79 


-0.41 






-0.49 


0.09 


3s 2 3p 2 P° /2 


3s 2 3d 2 D 3/2 


1260.42 


+0.39 


+0.32 


+0.49 


+0.40 


+0.38 






+0.38 


0.05 


3s 2 3p 2 P° /2 


3s 2 3d 2 D 3/2 


1265.02 


-0.28 


-0.39 


-0.58 


-0.68 


-0.36 






-0.40 


0.10 


3s 2 3d 2 D 5/2 


1264.73 


+0.61 


+0.70 


+0.61 


+0.57 


+0.63 






+0.62 


0.04 


Q„2n 2t-iO 

OS dp ±1/2 


o^o«2 2 t~> 

3s3p P1/2 


i i no o o 


+0.041 


-0.068 


+0.060 


+0.082 


+0.069 






+0.037 


0.047 


3s 3p P 3/2 


o _ o ,_2 2 td 

3s3p fi/2 


1 197.39 


-0.24 


-0.22 


-0.34 


-0.31 


-0.21 






-0.25 


0.04 


3s 2 3p 2 Pi/ 2 


3s3p 2 2 P 3/2 


1190.42 


-0.04 


-0.55 


-0.68 


-0.33 


-0.25 






-0.29 


0.10 


3s 2 3p 2 P° /2 


3s3p 2 2 P 3 / 2 


1194.50 


+0.39 


+0.52 


+0.39 


+0.38 


+0.48 






+0.28 


0.22 


3s3p 2 4 Pi/ 2 


3s 2 4p 2 P° /2 


o ct c\ c 


-6.24 


-6.45 


-7.39 


-6.90 


-9.34 






-6.56 


0.34 


3s3p 2 4 Pi/ 2 


3s 2 4p 2 P 3/2 


zdOI.od 


-6.49 


-7.33 


-6.54 


-6.03 


-7.99 






-6.50 


0.45 


3s3p 2 4 P 3/2 


3s 2 4p 2 P° /2 


2613.00 


-5 75 


-6 22 


-5 95 


-6 10 


-5.64 






-5.88 


0.21 


Ss^n 2 4 P, 
osop r 3/2 


Ss 2 4n 2 P° 


Q1 


-7.97 


-6.61 


-7.80 


-7.17 


-6.63 






-6.94 


0.40 


3s3p 2 4 P 5/2 


3s 2 4p 2 P° /2 


2620.90 


-4.87 


-5.56 


-5.00 


-5.33 


-4.63 






-4.97 


0.30 


3s3p 2 2 D 3/2 


3s 2 4p 2 P° /2 


3862.60 


-0.73 


-0.89 


-0.50 


-0.90 


-0.67 






-0.86 


0.28 


3s3p 2 2 D 3/2 


3s 2 4p 2 P° /2 


3853.66 


-1.43 


-1.58 


-1.20 


-1.60 


-1.37 






-1.55 


0.28 


3s3p 2 2 D 5/2 


3s 2 4p 2 P° /2 


3856.02 


-0.48 


-0.64 


-0.25 


-0.64 


-0.42 






-0.62 


0.32 


3s 2 4s 2 S 1/2 


3s 2 4p 2 P° /2 


6371.36 


-0.37 


-0.99 


-0.11 


-0.21 


-0.12 






-0.32 


0.20 


3s 2 4s 2 S 1/2 


3s 2 4p 2 P° /2 


6347.10 


-0.55 


-0.04 


-0.20 


-0.09 


+0.19 






-0.01 


0.19 
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data. Both, length and velocity (//-values, are taken with 
equal weights in the average. In the case of the recombi- 
nation transitions and the 1814A multiplet we also include 
in the averages the experimental data from Calamai et al. 
(1993) and Bergeson & Lawlcr (1993) respectively with 
twice as much weight as the theoretical values. It is noted 
that our calculation AST8 that uses the smaller configu- 
ration expansion is very well optimized on the lowest en- 
ergy terms, for instance the 3s 2 3p 2 P° and 3s3p 2 2 D terms, 
but it deteriorates rapidly for higher excitation terms. Fig. 
2 depicts our recommended ^/-values for every transition 
together with the statistical dispersion, which is probably 
representative of the true uncertainty. In Table 3 we show 
our results for transitions among the lowest 15 levels of 
the ion from the AST9, AST9, HFR, and MCDF calcu- 
lations. The AUTOSTRUCTURE results are given in length 
and velocity gauges and the MCDF results are given in 
the corresponding Coulomb and Babushkin gauges. In the 
last column of the table are our recommended values and 
estimated uncertainties 



3. Collision Strengths 

The collision strengths for electron impact excitation are 
computed with the ICFT Breit-Pauli i?-matrix package 
(bprm) based on the close-coupling approximation of Burke 
& Seaton (1971) whereby the wavefunctions for states of an 
A^-elcctron target and a colliding electron with total angu- 
lar momentum and parity Jtt are expanded in terms of the 
target eigenfunctions 

* 3 

The functions \% are vector coupled products of the target 
eigenfunctions and the angular components of the incident- 
electron functions, Fi(r) arc the radial part of the pro- 
jectile electron and A is an antisymmetrization operator. 
The functions $j are bound-type functions of the total sys- 
tem constructed with target orbitals; they are introduced 
to compensate for orthogonality conditions imposed on the 
Fi(r) and to improve short-range correlations. The Kohn 
variational principle gives rise to a set of coupled integro- 
diffcrential equations that are solved by i?-matrix tech- 
niques (Burke & Seaton 1971; Berrington et al. 1974, 1978, 
1987) within a box of radius, say, r < a. In the asymptotic 
region (r > a) exchange between the outer electron and the 
target ion can be neglected, and the wavefunctions can be 
approximated by Coulomb solutions. 

One-body Breit-Pauli relativistic corrections have been 
introduced in the i?-matrix suite by Scott & Burke (1980); 
Scott & Taylor (1982). Inter-channel coupling is equivalent 
to CI in the atomic structure context. 

Because of the large number of configurations and close 
coupling states in the representation of the target ion the 
scattering calculation had to be done in LS-coupling. Then, 
fine structure collision strengths were obtained by means 
of the intermediate coupling frame transformation (ICFT) 
method of Griffin et al. (1998). 

Since we have produced a large number of target ex- 
pansion to study the quality and the wavefunction it is in- 
teresting to see the effects of these various representations 
on the collision strengths. This comparison may be used to 




Fig. 5. Comparison of collision strengths in LS coupling 
for excitations from the ground term to the first five excite 
terms in Sin computed with approximations AST1 (left 
panels), AST8 (middle panels), and AST9 (right panels). 



asses the accuracy of the collision strengths. Thus, we per- 
formed three different scattering calculations using target 
orbitals from models AST1, AST8, and AST9. For the first 
two calculations we included only the lowest 12 LS terms 
in the closed coupling expansion, while for the calculation 
with the AST9 target, the most accurate and extensive, 
we included the lowest 43 LS terms. All calculations ex- 
plicitly include partial waves from states with L < 16 and 
multiplicity 1, 3, 5, and 7. The final collision strengths are 
produced with an energy resolution of 6 x 10~ 5 Ry. 

Figure 5 compares the LS-coupling collision strengths 
obtained from the target expansions AST1, AST8, and 
AST9 for excitation from the ground term 3s 2 3p 2 P° to 
the first four excited terms 3s3p 2 4 P, 3s3p 2 2 D, 3s 2 4s 2 S, 
and 3s3p 2 2 S. 

The results, in terms of Maxwcllian averaged effective 
collision strengths, of AST1, AST8, and AST9 are also com- 
pared in Table 4 for excitations from the ground term, 
3s 2 3p 2 P°, to the first five excited terms of Sill. Here, 
the third column presents the effective collision strengths 
from our best model AST9 and columns fourth and five 
present the percentage difference from these as obtained 
from models AST1 and AST8 respectively. These compar- 
isons are done for temperatures between 5000 K and 20,000 
K. Excitation rates for higher temperatures are of little 
practical interest because under these conditions the ionic 
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Table 4. Comparison of Maxwellian averaged collision 
strengths in LS-coupling for excitation from the ground 
term 3s 2 3p 2 P°. 



Upper term 


T (K) 


Model 9 


Difference (%) 








Model 1 


Model 8 


3s3p 2 4 P 


5000 


4.44E+0 


1.2 


5.7 




10000 


4.36E+0 


1.4 


5.9 




20000 


4.43E+0 


-2.4 


-5.2 


3s3p 2 2 D 


5000 


1.19E+1 


-8.2 


-39.8 




10000 


1.18E+1 


-3.1 


-36.4 




20000 


1.10E+1 


3.9 


-34.2 


3s 2 4s 2 S 


5000 


2.17E+0 


11.0 


25.3 




10000 


2.18E+0 


16.3 


27.0 




20000 


2.21E+0 


15.2 


33.1 


3s3 P 2 2 S 


5000 


3.32E+0 


-38.1 


-52.4 




10000 


3.08E+0 


-32.7 


-54.7 




20000 


2.95E+0 


-24.4 


-59.2 



fraction in Sin is too small. The are two general charac- 
teristics of collision strengths that could lead to variations 
in their thermal averages: (i) the background and reso- 
nances in the near threshold region, which determine the 
low temperature (AE/kT < 1) Maxwellian averaged col- 
lision strengths and depend on the coupling of the target 
ion with the continuum; and (ii) the slope of the continuum 
towards high energies, which in dipole allowed transitions 
depend linearly on the oscillator strength and consequently 
on the quality of the target representation. It is interest- 
ing to see that the results from AST1 and AST9 agree 
rather well in terms of the qualitative shape of the collision 
strengths and the quantitative Maxwellian averaged rates, 
within ~ 10%, except for the 3s 2 3p 2 P° - 3s3p 2 2 S transi- 
tion. This is in contrast to the large differences in gf-values 
for essentially all transitions among levels of the TV-electron 
Si II target. The good agreement in collision strengths is be- 
cause these are dominated by bound-continuum couplings 
of the (N + l)-electron system, which seems to be well rep- 
resented by the large amount CI in the close coupling ex- 
pansion 

On the other hand, the target expansion AST8 does 
yield a different lay out of resonances on the collision 
strengths. In addition, AST8 yields steeper rises in the 
background collision strengths for excitation to the 3s 2 4s 
2 S level and higher levels. This translates into effective col- 
lision strengths from AST8 that are systematically lower 
than those from AST9 by as much as ~ 50%. 

It is expected that our results from the AST9 expansion 
should be the most accurate due to the good quality of the 
target representation and the large amount of CI included 
for the (N+l)-electron system. Nonetheless, the comparison 
with results of various target representation allow us to 
asses the uncertainties in the collision strengths to ~ 10%. 

Fig. 6 shows the collision strengths for a sample of for- 
bidden, intercombination, and dipole transitions. We chose 
the same transitions as in Tayal (2008) for direct compar- 
ison with that work. Although the collision strengths from 
both calculations are roughly similar, there are quantitative 
differences between the two sets in terms of resonance struc- 
tures and absolute level of the background. It is apparent 
that the present collision strengths should yield Maxwellian 
averages somewhat lower than those of Tayal. 



12 p — 1 — 1 — 1 — i — n 1 — 1 — i — 1 — 1 — 1 — r 




0.5 0.6 0.7 0.8 0.9 1 



Electron energy (Ry) 

Fig. 6. Collision strengths in JJ-coupling for transitions 
3s 2 3p 2 P° /2 - 3s 2 3p 2 P° /2 (top panel), 3s 2 3p 2 P° /2 - 

3s3p 2 4 P 3/2 (middle panel), and 3s 2 3p 2 P° /2 - 3s3p 2 2 D 5/2 
(bottomr panel). 

In Table 5 we compare the present fine structure 
Maxwellian averaged collision strengths with those of Tayal 
(2008) and Dufton & Kingston (1991). For most transitions 
our results are ~ 30% lower than those of Tayal, while the 
results of Dufton & Kingston lie inbetween those two. That 
our results are somewhat lower than those of Dufton & 
Kingston can be understood from the much larger close 
coupling expansion used in the present work. As more scat- 
tering channels are open in the calculation, the flow of 
electrons is redistributed and the collision strength among 
low lying levels tends to converge to lower values. The 
source of the differences found with respect to Tayal are 
less clear, and seems to be due to differences in computa- 
tional approach used, i.e. between the orthogonal and non- 
orthogonal R-matrix methods. 

4. Conclusions 

We have carried out extensive calculations of transitions 
rates and collision strengths for electron impact excitation 
for the lowest 12 levels of the astrophysically important Si n 
ion. 

In the calculation of radiative data, we paid special at- 
tention to the weak dipole allowed and intercombination 
transitions, which are of particular interest for plasma di- 
agnostics. Determination of accurate data for these transi- 
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Table 5. Comparison of Maxwellian averaged collision 
strengths in JJ-coupling from present calculation (Present), 
Tayal (2008; Tayal), and Dufton & Kingston (1991; DK). 



Uppci 


■ level 


T (K) 


Present 


Tayal 


DK 


3s^3p 


2-po 
r 3/2 


5000 


4.55 


6.19 


5.60 




10000 


4.45 


6.09 


5.70 






20000 


4.42 


5.97 


5.77 


3s3p 2 


4 Pl/2 


5000 


0.401 


0.512 


0.550 




10000 


0.398 


0.515 


0.516 






20000 


0.392 


0.502 


0.466 


3s3p 2 


4 P 3 /2 


5000 


0.612 


0.812 


0.832 




10000 


0.609 


0.789 


0.780 






20000 


0.602 


0.769 


0.706 


3s3p 2 


4 P 5 /2 


5000 


0.441 


0.615 


0.571 




10000 


0.458 


0.595 


0.534 






20000 


0.477 


0.589 


0.488 


3s3p 2 


2 D 1/2 


5000 


1.82 


2.77 


2.76 




10000 


1.82 


2.74 


2.74 






20000 


1.75 


2.50 


2.58 


3s3p 2 


2 P>3/2 


5000 


2.05 


2.94 


2.45 




10000 


2.14 


2.98 


2.44 






20000 


2.14 


2.80 


2.30 


3s 2 4s 


2 Sl/2 


5000 


0.910 


1.02 


1.24 




10000 


0.865 


1.06 


1.20 






20000 


0.857 


0.979 


1.04 


3s3p 2 


2 Sl/2 


5000 


0.887 


0.102 


0.716 




10000 


0.899 


0.988 


0.840 






20000 


0.916 


0.988 


0.902 



tions is particularly challenging, therefore we studied the 
effects of valence- valence, valence-core, and core-core inter- 
actions with three different methods, i.e. MCDF, HFR, and 
the central potential method implemented in AUTOSTRUC- 
ture. With MCDF we could only include valence-valence 
correlations, as opening of the n=2 core resulted in a large 
number of states that could not be managed with the com- 
puter code GRASP. For this reason, we were able to obtain 
accurate transition rates for the 3s 2 3p 2 P° - 3s3p 2 4 P inter- 
combination transitions only. Both HFR and AUTOSTRUC- 
ture allowed us to investigate valence-core and core-core 
interaction by building very large configuration expansions. 
The accuracy of AUTOSTRUCTURE calculations were signifi- 
cantly improved by the use of the modified TFDA potential 
of Bautista (2008) and a new technique for optimizing the 
variational parameters of this potential. This optimization 
technique takes into account the differences between length 
and velocity gauges of the g/-values, in addition to the 
accuracy of predicted energy levels. Our most accurate gf- 
values were then compared with previous theoretical and 
experimental determinations. From these comparisons we 
derive a recommended set of (//-values and estimate their 
uncertainties. 

We then proceed to compute electron impact excitation 
collision strengths with the R-matrix method. We do so by 
using various of the targets representations made from the 
previous calculations and compare the results. This allows 
us to identify the physical effects that affect the accuracy 
of the computed collision strengths. We also compare our 
results with those of previous calculations. The present re- 
sults agree reasonably well with those of Dufton & Kingston 
(1991), who also used an orthogonal R-matrix method but 



with a much smaller close coupling expansion. On the other 
hand, the present results for Maxwellian averaged colli- 
sion strengths are systematically lower than those of the 
recent calculation of Tayal (2008) using a non-orthogonal 
R-matrix method. The reasons for these differences, that 
typically amount to ~ 30%, are unclear. We argue that 
the source of this difference could be in the non-orthogonal 
R-matrix approach. 

Tables 6 and 7 containing the present gf-values, A- 
values, and effective collision strengths are available in elec- 
tronic form. The data and atomic model are also to become 
available through the TIPTOP 1 database and the XSTAR 
database Bautista & Kallman (2001). 
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